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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  the  establishment  of  receiving,  processing,  and  interpretation  ca¬ 
pabilities  for  spaceborne  along-track  interferometric  synthetic  aperture  radar  (along-track  InSAR,  ATI) 
data  at  the  University  of  Miami's  Center  for  Southeastern  Tropical  Advanced  Remote  Sensing 
(CSTARS)  for  high-resolution  surface  current  field  mapping  with  special  emphasis  on  applications  in 
estuarine  environments  and  rivers. 

OBJECTIVES 

Within  the  project  period  of  24  months,  which  started  officially  in  March  2009,  we  intended  to 

►  acquire  series  of  ATI  images  of  selected  test  sites; 

►  implement  ATI  raw  data  processing  capabilities  at  CSTARS; 

►  develop  a  robust  current  retrieval  algorithm  for  rivers,  which  corrects  the  ATI-derived  Doppler  ve¬ 
locities  for  contributions  of  wave  motions; 

►  compare  the  results  with  available  reference  data; 

►  evaluate  the  data  quality  in  view  of  known  user  requirements;  and 

►  evaluate  the  potential  of  dedicated  ATI  satellite  missions  with  optimized  instrument,  platform,  and 
orbit  parameters  for  nearshore,  estuarine,  and  river  observations. 

The  main  satellite  to  be  used  is  the  German  TerraS AR-X,  but  we  were  optimistic  to  obtain  additional 
ATI  data  from  other  satellites,  such  as  TanDEM-X,  RADARSAT-2,  and  the  Italian  COSMO-SkyMed 
satellite  constellation. 

APPROACH 

The  project  consists  of  five  main  tasks,  which  are  defined  as  follows. 

Task  1:  Data  Acquisition.  Test  sites  of  interest  were  to  be  selected,  taking  into  account  requirements  of 
ONR  and  its  existing  programs,  the  suitability  of  the  test  sites  for  ATI  data  acquisitions  and  for  the 
tests  to  be  performed,  and  the  availability  of  reference  data  for  validation  of  the  results.  TerraSAR-X 
ATI  images  of  the  test  sites  were  to  be  acquired  repeatedly  during  a  period  of  several  months.  At  the 
beginning  this  was  expected  to  be  done  in  close  collaboration  with  Hartmut  Runge  and  Ulrich  Stein- 
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brecher  at  DLR  Oberpfaffenhofen  (Germany),  who  would  order  the  experimental  data  products  and 
program  the  instrument  for  the  individual  data  acquisitions,  respectively.  At  a  later  stage,  the  ordering 
of  ATI  data  acquisitions  with  TerraSAR-X  would  be  possible  through  the  Web  interface  that  is  also 
used  for  ordering  conventional  data  products  by  regular  Pis  of  TerraSAR-X.  Furthermore,  CSTARS 
would  build  up  its  own  data  ordering  system  for  TerraSAR-X  and  obtain  direct  access  to  the  satellite  in 
2010.  As  a  secondary  activity,  Task  1  includes  the  collection  of  reference  data  from  in-situ  measure¬ 
ments,  numerical  models,  and  existing  data  bases.  It  was  planned  to  obtain  such  data  from  colleagues 
at  other  research  institutions  and  agencies. 

Task  2:  Raw  Data  Processing.  The  ATI  raw  data  need  to  be  processed  into  complex  SAR  images,  and 
some  additional  processing  steps  (co-registration  etc.)  need  to  be  applied  to  obtain  the  phase  images 
required  for  current  retrievals.  While  this  work  was  originally  done  by  Steffen  Suchandt  at  DLR  Ob¬ 
erpfaffenhofen  (Germany),  we  were  planning  to  obtain  his  software  for  use  at  CSTARS  in  the  course 
of  this  project.  The  project  was  expected  to  make  a  key  contribution  to  the  establishment  of  ATI  data 
receiving  and  processing  capabilities  at  CSTARS. 

Task  3:  Current  Retrieval.  The  conversion  of  ATI  phase  images  into  actual  line-of-sight  current  fields 
requires  an  application  of  advanced  filtering  techniques  and  a  correction  of  detected  Doppler  velocities 
for  contributions  of  wave  motions  and  their  spatial  variations  due  to  wave-current  interaction,  as  theo¬ 
retically  discussed  and  partly  demonstrated  in  previous  publications  [1][2][3],  The  ATI-derived  current 
fields  from  the  Shuttle  Radar  Topography  Mission  (SRTM)  in  [2]  and  [3]  were  not  corrected  by  this 
method  because  they  were  not  absolutely  calibrated,  thus  the  mean  contribution  of  wave  motions  could 
not  be  separated  from  other  contributions  that  were  not  known  theoretically.  The  absolute  bias  of  the 
detected  velocities  had  to  be  tuned  by  different  methods.  This  is  not  necessary  with  TerraSAR-X  data, 
but  some  theoretical  model  and  algorithm  improvements  are  required  for  river  applications,  since  the 
computation  of  spatially  varying  wave  contributions  to  Doppler  velocities  is  a  different  problem  for  a 
river  than  for  open  waters.  The  model  must  be  upgraded  to  account  for  shallow-water  effects  and 
boundary  effects  such  as  wave  refraction,  reflection,  and  dissipation  at  river  banks.  All  filtering  and 
correction  of  the  ATI  data  and  the  development  of  model  improvements  was  to  be  done  by  PI  Roland 
Romeiser.  In  addition,  an  algorithm  for  direct  wind  retrievals  from  the  amplitudes  of  TerraSAR-X 
SAR  /  ATI  images  was  to  be  implemented  to  eliminate  the  need  for  wind  data  from  external  sources. 
This  was  to  be  done  in  collaboration  with  Jochen  Horstmann  of  NURC,  La  Spezia  (Italy). 

Task  4:  Evaluation.  The  ATI-derived  surface  current  fields  must  be  checked  for  consistency  and  com¬ 
pared  with  available  reference  data  to  determine  the  achievable  data  quality  under  varying  conditions 
(wind  speed  /  surface  roughness,  river  width,  angles  between  flow  direction,  wind  direction,  and  radar 
look  direction,  etc.).  Based  on  the  results,  the  suitability  of  the  method  for  various  applications  will  be 
discussed,  and  the  potential  of  dedicated  spacebome  ATI  missions  with  optimized  instrument  and  orbit 
parameters  for  ocean  and  river  applications  will  be  evaluated. 

Task  5:  Dissemination.  The  results  obtained  at  different  stages  of  the  project  were  to  be  disseminated 
through  a  project  website,  periodic  progress  reports  to  ONR  (such  as  this  one),  presentations  at  infor¬ 
mal  meetings  as  well  as  national  and  international  conferences,  and  at  least  two  publications  in  peer- 
reviewed  journals.  To  ensure  an  efficient  collaboration  and  data  exchange  with  existing  ONR  program 
teams  and  other  potential  partners  and  users,  we  were  planning  to  contact  these  colleagues  in  an  early 
stage  of  the  project  and  inform  them  about  our  objectives,  plans,  and  needs. 
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WORK  COMPLETED 


Task  1:  Data  Acquisition.  First  TerraS  AR-X  ATI  images  of  ocean  and  river  scenes  were  acquired  in 
Aperture  Switching  (AS)  mode  in  spring  and  summer  2008,  including  six  images  of  the  Elbe  river 
(Germany),  which  were  made  available  to  us  in  December  2008.  In  spring  2009,  another  set  of  images 
of  the  Elbe  river  was  acquired  in  Dual  Receive  Antenna  (DRA)  mode.  In  Fall  2009,  the  acquisition  of 
AS-mode  images  was  continued.  The  Elbe  river  was  selected  as  first  test  site  a  long  time  ago  because 
we  had  experiences  with  it  from  the  analysis  of  SRTM  data  [3]  and  it  was  easy  to  obtain  additional  ref¬ 
erence  data  from  governmental  agencies.  After  the  Elbe  river,  the  Lena  river  in  Siberia  was  selected  as 
a  second  river  for  routine  data  acquisitions  during  a  longer  period.  The  Lena  river  has  strong  seasonal 
variabilities  in  its  runoff  (minimum  1400  m3/s  in  April,  maximum  74,000  m3/s  in  June),  and  it  is  cov¬ 
ered  with  ice  for  several  months  in  winter.  Reference  data  are  available  from  a  station.  In  April  and 
May  2010,  a  series  of  DRA-mode  images  was  acquired  over  Skagit  Bay  (WA),  Columbia  River  (OR  / 
WA),  St.  Johns  River  (FL),  Amazon  (Brazil),  and  a  few  open-ocean  test  sites  and  other  test  sites.  The 
processing  of  these  data  is  under  way.  Receiving  capabilities  for  TerraS  AR-X,  COSMO-SkyMed,  and 
RADARSAT-2  at  CSTARS  were  implemented  as  planned,  but  the  existing  generation  of  COSMO- 
SkyMed  satellites  does  not  have  the  ATI  data  acquisition  capabilities  that  were  discussed  in  some  con¬ 
cept  studies,  thus  we  have  not  obtained  ATI  data  from  other  satellites  than  TerraS  AR-X.  A  proposal 
for  some  special  ATI  and  combined  ATI-XTI  (along-track  and  cross-track  interferometry)  experiments 
with  the  newly  available  TerraS  AR-X  /  TanDEM-X  constellation  (the  TanDEM-X  satellite  was 
launched  as  a  companion  to  TerraS  AR-X  in  summer  2010)  was  submitted  to  DLR  and  is  currently  un¬ 
der  review. 

Task  2:  Raw  Data  Processing.  The  first  six  TerraS  AR-X  AS-mode  images  of  the  Elbe  river  were 
processed  in  December  2008.  An  analysis  revealed  some  artifacts  in  the  data  that  could  be  attributed  to 
suboptimal  processing  parameters.  Improved  results  were  obtained  in  September  2009.  The  original 
processing  of  DRA-mode  images  from  spring  2009  did  not  lead  to  satisfactory  results,  but  improved 
procedures  permitted  current  retrievals  of  quite  high  quality  in  spring  2010.  Examples  of  results  ob¬ 
tained  in  2010  are  shown  at  the  end  of  this  document. 

Task  3:  Current  Retrieval.  Using  the  first  six  TerraSAR-X  AS-mode  images  of  the  Elbe  river,  specific 
filtering  and  correction  algorithms  for  TerraSAR-X  data  were  developed.  These  algorithms  can  identi¬ 
fy  and  eliminate  phase  signatures  of  ships  and  other  strong  scatterers  in  the  water  with  nonrepresenta¬ 
tive  velocities,  reduce  the  effect  of  azimuth  ambiguities  (ghost  images  of  strong  scatterers  on  land  over 
water),  correct  the  absolute  phase  calibration  on  the  basis  of  phases  measured  over  land,  reduce  re¬ 
maining  phase  noise  by  filtering,  convert,  ATI  phases  into  radial  velocities,  and  correct  the  velocities 
for  contributions  of  wave  motions.  Almost  the  same  algorithms  can  be  applied  to  the  DRA-mode  im¬ 
ages  from  spring  2009  and  spring  2010.  Some  specific  ATI  model  modifications  for  river  applications 
have  been  implemented,  but  they  need  to  be  tested  with  data  from  a  variety  of  different  test  sites  and 
current  and  wind  conditions.  This  is  an  ongoing  task. 

Task  4:  Evaluation.  Aside  from  some  uncertainties  about  the  best  choice  of  model  settings  for  the  cor¬ 
rection  of  ATI  data  for  contributions  of  wave  motions,  our  latest  results  show  good  agreement  between 
ATI-derived  currents  and  reference  currents  and  a  quite  consistent  data  quality  for  all  images  (and  for 
different  regions  within  each  image).  A  final  evaluation  will  be  presented  when  all  data  from  the  DRA- 
mode  experiment  in  spring  2010  have  been  processed  and  the  results  have  been  analyzed. 
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Task  5:  Dissemination.  The  first  TerraS AR-X  results  obtained  in  2009  were  presented  in  a  press  re¬ 
lease  of  RSMAS  and  reproduced  on  several  web  sites.  A  full  peer-reviewed  paper  was  accepted  for 
publication  in  IEEE  Trans.  Geoscience  Remote  Sensing  in  August  2009  and  published  in  February 
2010  [4],  Furthermore,  the  early  results  were  presented  at  Oceans  2009  in  Bremen,  Germany  [5],  at  the 
2009  ONR  Physical  Oceanography  Review  Symposium  in  Chicago,  and  at  the  2009  International 
Geoscience  and  Remote  Sensing  Symposium  (IGARSS)  in  Cape  Town,  South  Africa  [6],  More  results 
were  presented  at  Oceans  From  Space  in  Venice,  Italy  [7],  The  Living  Planet  Symposium  of  the  Euro¬ 
pean  Space  Agency  in  Bergen,  Norway  [8],  and  IGARSS  2010  in  Honolulu  [9], 

RESULTS 

Some  current  fields  (and  ship  speed  results)  obtained  for  Elbe  river  images  are  shown  in  Figs.  1-4.  The 
data  quality  is  good  and  consistent  with  our  expectations.  More  such  results  from  the  2010  DRA-mode 
campaign  will  become  available  soon.  We  know  now  how  to  process  and  correct  such  data,  and  we 
have  receiving  capabilities  for  TerraSAR-X  at  CSTARS.  We  are  planning  to  do  current  measurements 
by  TerraSAR-X  over  some  more  rivers  and  coastal  areas  in  a  follow-on  project,  in  which  we  would 
like  to  develop  algorithms  for  deriving  higher-level  data  products  such  as  river  runoff  estimates  and 
bathymetric  maps  of  coastal  waters. 

IMPACT/APPLICATIONS 

Techniques  for  satellite-based  current  measurements  in  rivers  and  coastal  areas  are  valuable  for  a  va¬ 
riety  of  applications  in  oceanography,  climate  research,  hydraulic  engineering,  drinking  water  and 
flood  risk  management,  and  other  disciplines.  For  DoD  operations,  the  navigability  of  rivers  and  the 
monitoring  and  prediction  of  flooding  events  and  changes  in  the  water  supply  for  major  cities  and  agri¬ 
cultural  regions  are  of  interest.  The  use  of  satellite-based  monitoring  techniques  is  the  most  promising 
approach. 

RELATED  PROJECTS 

The  techniques  developed  here  will  be  used  for  other  ONR-funded  projects;  such  as  the  DRI  "Inlet  and 
River  Mouth  Dynamics"  in  which  we  are  involved.  Plans  for  more  ATI-based  projects  at  RSMAS  / 
CSTARS  exist. 
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Fig.  1:  Line-of-sight  surface  current  fields  in  the  Elbe  river  as  obtained  for  the  six  first  test  images 
from  spring  2008  with  the  latest  processing  software.  White  lines  indicate  boundaries  of  the  main 
test  area  for  which  statistical  analyses  have  been  carried  out.  Signatures  over  land  show  unfiltered 
zero-mean  residual  Doppler  velocities.  As  a  result  of filtering,  the  current  fields  in  the  river  look 
very  smooth.  Mean  currents  in  the  six  cases  range  from  about  —1.9  to  1.3  m/s.  Artifacts  that  were 
clearly  visible  in  the  first  published  versions  of  these  images  could  be  eliminated  by  optimizing  filter 

parameters  in  the  SAR  raw  data  processing  algorithm. 
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Fig.  2:  Line-of-sight  surface  current  fields  in  the  Elbe  river  according  to  numerical  model  UnTRIM 
of  the  German  Federal  Waterways  Engineering  and  Research  Institute  (RAW),  for  six  tidal  phases 
corresponding  to  the  tidal  phases  at  the  times  of  the  TerraSAR-X  overpasses.  The  current  fields  look 

similar  to  the  ones  obtained  from  TerraSAR-X. 
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Fig.  3:  Diagram  showing  temporal  variations  of  spatial  mean  line-of-sight  currents  in  the  mean  test 
area  as  function  of  tidal  phase,  according  to  UnTRIM  (gray  bars)  and  TerraSAR-X.  Blue  diamonds 
show  uncorrected  Doppler  velocities;  green  triangles  show  Doppler  velocities  after  phase 
recalibration;  red  lines  between  triangles  and  red  squares  indicate  theoretical  maximum  of 
corrections  for  wave  contributions  according  to  our  model,  and  red  asterisks  show  best  estimates  of 
mean  currents  after  applying  50%  of  these  corrections.  The  final  results  are  in  pretty  good 

agreement  with  UnTRIM. 
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Fig.  4:  Two  high-resolution  SAR  intensity  images  (top)  and  derived  horizontal  line-of-sight 
velocities  (bottom)  for  a  different  part  of  the  Elbe  river  where  DRA-mode  images  were  acquired  in 
spring  2009.  The  velocity  images  are  similar  to  the  ones  of  Fig.  1,  but  the  data  quality  is  clearly 
better  because  of  a  higher  original  pixel  resolution  and  a  better  overall  data  quality  of  DRA-mode 
data.  These  images  were  processed  such  that  ships  (bright  spots  in  the  intensity  images)  were 
included  in  the  ATI  data  analysis;  their  speeds  are  shown  with  a  different  color  coding.  While  the 
mean  line-of-sight  currents  in  these  images  are  on  the  order  of  1  to  2  m/s  to  the  left,  the  line-of-sight 
speeds  of  9  ships  in  the  left  and  7  ships  in  the  right  image  range  from  about  7  m/s  to  the  left  to  7  m/s 

to  the  right. 
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